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ABSTRACT 

Aims. We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, 
VETS 450 (09.7H-Ib + 07::) and VETS 652 (B1 Ib + 09: HI:). 

Methods. We perform model-atmosphere analyses to characterise the photospheric properties of both members of each 
binary (denoting the ‘primary’ as the spectroscopically more conspicuous component). Radial velocities and optical 
photometry are used to estimate the binary-system parameters. 

Results. We estimate Tes = 27 kK, \ogg = 2.9 (cgs) for the VETS 450 primary spectrum (34 kK, 3.6: for the secondary 
spectrum); and Tes = 22 kK, loggr = 2.8 for the VETS 652 primary spectrum (35 kK, 3.7: for the secondary spectrum). 
Both primaries show surface nitrogen enrichments (of more than 1 dex for VETS 652), and probable moderate oxygen 
depletions relative to reference LMC abundances. We determine orbital periods of 6.89 d and 8.59 d for VETS 450 and 
VETS 652, respectively, and argue that the primaries must be close to Hlling their Roche lobes. Supposing this to be 
the case, we estimate component masses in the range ~20-50M©. 

Conclusions. The secondary spectra are associated with the more massive components, suggesting that both systems 
are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the 
spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra 
are contaminated by (or arise in) accretion disks. 

Key words, stars: early-type - binaries: spectroscopic - stars: variable: general - stars: fundamental parameters - stars: 
individual (VETS 450, VETS 652) 


1. Introduction 

Massive, luminous stars are of interest for the role that they 
play in galactic chemical evolution; the environmental im¬ 
pact they have through mechanical and radiative energy 
input to their surroundings; and as tracers of recent star 
formation. However, while there have been considerable ad¬ 
vances in modelling their spectra, direct determinations of 
their fundamental parameters are relatively few, because of 
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* Based on observations obtained at the European Southern 
Observatory Very Large Telescope (VLT) as part of programmes 


182.D-0222 and 090.D-0323 


the scarcity o f suitable doub le-lined eclipsing binary sys¬ 
tems (cf., e.g., Bonanos||2009 and references therein). 

Multi-epoch spectros copy from the VLT-FLAMES 
Tarantula Survey (VETS; Evans et al.|201I ) of the OB-star 
population of 30 Doradus, in the Large Magellanic Cloud 
(LMC), has led to the discovery of a number of systems 
showing radial-velocity va riations that appear to be consis¬ 
tent with binary motion (Sana et al. 2013 Dunstall et al. 
20151. These systems offer an important opportunity to bet¬ 


ter understand the physical properties of stars in the upper 
Hertzsprung-Rus sell diagram, as exem plified by the VETS 
study of R139 by Taylor et al. (20111. 


Here we discuss two newly identihed double-lined radial- 
velocity variables discov ered in the V ETS: nos. 450 and 
652 (stars 50 and 5 in Melnick 1985). The primary ve- 
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Fig. 1: The central region of 30 Dor, showing R136 and the 
locations of VFTS 450 and VFTS 652. 


locity amplitudes are among the largest measured in the 
VFTS dataset, and the number of spectroscopic observa¬ 
tions available makes it possible to undertake full orbital 
analyses (rather than mere detections of variability). Each 
system also shows orbital photometric variability. Both are 
relatively close to the core of 30 Dor, with radial distances 
of 0147 and 0183 from R136 (Fig. [^, corresponding to pro¬ 
jected distances of 6.8 and 12.0 pc, respectively. 

The paper is organized as follows: Section summa¬ 
rizes the data, and the binary characteristics are examined 
in Section A model-atmosphere analysis is described in 
Section]^ and simple models of the systems are constructed 
in Section]^ Throughout the paper we adopt the conven¬ 
tion that the ‘primary’ in each system is the star with the 
stronger optical absorption-line spectrum (though we shall 
argue that this is probably not the more massive compo¬ 
nent). 


Table 1: Characteristics of spectroscopic observations. 


Wavelength 

setting 

A range 

(A) 

Resolving 
power R 

Typical 

S:N 

LR02 

3960-4560 

7000 

100-300 

LR03 

4505-5050 

8 500 

150 

UVES 520 

4175-5155, 

5240-6200 

53 000 

20-40 


2. Observations 

2.1. Optical Spectroscopy 

Initial spectroscopic data were obtained as part of the 
VFTS (Evans et al. 2011), using the Fi bre Large Array 
Mult i-Element Spectrograph (FLAMES; Pasquini et al. 
20021 on the Very Large Telescope, primarily with the 
Giraffe spectrograph, but with supplementary data from 
the Ultraviolet and Visual Echelle Spectrograph (UVES). 


These observations were obtained in the 2008/9 and 
2009/10 observing seasons; additional Giraffe spectroscopy 
was secured as part of a binary-monitoring campaign of 
VFTS targets between 2012 Oct and 2013 Mar. 

Table □ summarizes the basic instrumental characteris¬ 
tics; a full account of the obs ervations and data reduction 
is given by Evans et’ah] ( |2011[ ). Logs of the individual blue- 
region spectra, which are the principal focus of this paper, 
are given in Tables [^and[^ (on-line). Representative spectra 
are shown in Fig. 


2.1.1. Spectral types 

Spectral types previously determined from the VFTS spec¬ 
tra are 09.7 III:-b 07:: and B2 Ip -b09 III: (VFT S 450, 


VFTS 652; Walborn et al. 2014). Melnick 


09.5 I and 09.5 V pec (“binary?”) for Vl'T S 


1985) gives 
50 and 652, 

respectively, while Walborn & Blades (1997) report ON9:1 
and B2 Ib. 

Our review of the more extensive dataset discussed here, 
including examination of the disentangled component spec- 
tra presente d in Section broadly supports the Walborn 
et al. (2014) classihcations, but the clear presence of Siiv 
A4089 and A4116 in the primary spectrum of VFTS 652 
leads us to revise its classification to B1 Ib. The crucial 
He I A447I classihcation line suffers significant nebular con¬ 
tamination in the disentangled secondary spectra, admit¬ 
ting the possibility of an 08 (or, conceivably, 07) secondary 
spectrum for this target. 

Our rectihcation of the VFTS 450 spectra leaves a 
broad, shallow emission feature spanning AA4640, 4686 
(C iii/N III, He II; Fig. [2]). We have investigated, and re¬ 
jected, possible instrumental origins, including contamina¬ 
tion by the nearby WR star Brey 79 (3'.'5 distant). While 
such features are not widely reported, and are easily over¬ 
looked, they are not unprecedent ed in late-0 supergiants 
(e.g., a Cam, 09 la; Wilson|[l958 ); this suggests the possi¬ 
bility of a brighter luminosity class for the VFTS 450 pri¬ 
mary than previously inferred from VFTS data. The inten¬ 
sity of the Si IV lines compared to He i A4026 als o indicates 
a som ewhat more luminous type (cf. Table 6 of Sota et al. 
2011). The arbitrary intensity scaling of the disentangled 


spectrum hampers a precise assignment, but we revise the 
previous classification for the VFTS 450 primary spectrum 
to 09.7 II-IbFlOur adopted spectral types are incorporated 
into Tables^ 


2.1.2. Ha spectra 

We have Ha observations, shown in Fig. at a single 
epoch for each system. These spectra suffer from strong 
nebular contamination which is poorly corrected by stan¬ 
dard sky subtraction, but nevertheless each star clearly 
shows broad, double-peaked intrinsic emission. Although 
the single-epoch spectra may not be representative of typ¬ 
ical behaviour, this emission morphology is characteristic 
of interacting binaries, rather than typical OB-star stellar- 
wind P-Cygni prohles. Peak-to-peak separations are ~560 
and 420 km s“^ for VFTS 450 and 652, respectively, with 
full widths at continuum level about twice those values. 

^ We recall the convention that ‘Il-Ib’ is to be read as indi¬ 
cating a range of uncertainty, whereas ‘Ib-II’ would indicate a 
precise interpolated luminosity class. 
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Table 4: Basic observed properties. 


VFTS 

Spectral Type 

V 

to 

1 

J 

H 

Ks 

[3.6] 

[4.5] 

[5.8] 

< /c > 

<V-Ic> 

450 

09.7II-IbT07:: 

13.60 

0.20 

13.08 

12.91 

12.89 

11.00 

10.65 

10.52 

13.26 

T0.19 

652 

B1 IbT09: III: 

13.88 

0.20 

13.40 

13.28 

13.22 

12.97 


12.73 

13.63 

T0.36 


Notes. Photometry follows |Bonanos et al.|(2009|). The primary source for the B, V photometry is ISelman et al. (|1999|); their 
observations were obtained over only ~10 ininutes, so the {B — V) colours are insensitive to orbital variability. JHK results are 
from IRSF (the InfraRed Survey Facility; |Kato et al.||2007|), and the mid-IR photometry from Spitzer ‘SAGE’ Legacy Science 
Program ( Meixner et al.|2006 |. Measurement uncertainties are <0™05, excepting the [5.8] magnitude for VFTS 450 (TOT 15), but 
both stars are variable with amplitudes of ~0T2 (Fig. |5|. The last two columns are average results from our OGLE photometry 

(pnj). 



Wavelength (A) 


Fig. 2: Rectified blue-region spectra of VFTS 450 and 652, velocity shifted to the rest frame of the primary. The data 
are LR02 and LR03 spectra from MJD 54748, 54810 (VFTS 450; ~ 0.12 from the circular-orbit ephemeris in Table 

and 56294, 54808 (VFTS 652, cj) ~ 0.66), merged at ~4560A. Secondary spectra are offset by ca. T260, —150 km s“^ 
(VFTS 450, 652, respectively). Narrow Balmer emission is nebular. 


2.2. Photometry 
2.2.1. OGLE photometry 


We have OGLE /c-b and photometry from pha ses III and 
IV of the project (cf. [Udalski et al. 2008 20151, spanning 
2001 October to 2009 April and 2010 March to 2014 March, 
respectively. Each dataset for each star consists of ~400 
observations. We have also examined the sparser OGLE-III 
V-band data. 

VFTS 450 is located in a high-background region, 
and, as noted above, is only ^375 arcsec from Brey 79 
{Ic — 12.7). The standard OGLE-III Difference Image 
Analysis (DIA) pipeline is not optimal under these circum¬ 
stances, and w e found that a profile- fitting extraction, us¬ 
ing DoPhot ( Schechter et al.| 1993[ ), resulted in reduced 
scatter. Furthermore, the OGLE-IV Ic photometry for this 
target is ^0T17 brighter than the OGLE-III data (regard¬ 
less of extraction method). This is probably a consequence 
of the high background; while it is difficult to be certain, we 


believe the OGLE-IV normalization to be the more reliable. 
Neither issue arises in the VFTS 652 results. 


Both stars show orbital photometric variability, with full 
amplitudes of ^0T2. Periods were determined by using a 


date-compensated discrete fourier transform (Ferraz-Mello 


19811, augmented with least-squares fitting of a double sine 


wave. Results are included in Table There are no signif¬ 
icant differences in periods determined from the OGLE- 
III, OGLE-IV, and combined datasets (Table gives re¬ 
sults from the combined /c-band data). The phased OGLE 
photometry is shown in Fig. The rms scatter about the 
mean curve for VFTS 450, ~0T023, is consistent with the 
probable measurement uncertainties, but the larger scatter 
for VFTS 652, suggests significant intrinsic vari¬ 

ability. 


2.2.2. IR photometry 


Representative visual-IR mag nitudes for both s tars, 
adapted from the compilation by Bonanos et al. (2009), are 
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Wavelength (A) 


Fig. 3: Ha spectra, labelled with MJDs of mid-observation. 
The subtracted sky spectrum is shown for reference, and il¬ 
lustrates the nebular contamination, which varies on small 
spatial scales; correction for this nebular emission in the 
stellar spectra is generally poor, in particular in the resid¬ 
ual core Ha emission, although the extended double-peaked 
emission is real. 



J - [3.6] 


Fig. 4: Colour-magnitude diagram, comparing VFTS 450 
and 652 to selec ted LMC emission-li ne stars and OBA 
supergiants (after Bonanos et al.|[2009). 


listed in Tableand are included in a J — [3.6], [3.6] colour- 
magnitude diagram of luminous LMC sources in Fig. [^ 
VFTS 652 lies at the red edge of the distribution of normal 
OBA supergiants in this figure (though this displacement 
from the main grouping could possibly arise from differ¬ 
ent phase sampling at J and [3.6] of the orbital photo¬ 
metric variability discussed in Section 3.21. However, at 
these wavelengths VFTS 450 has a substantial IR excess, 
intermediate between those of typical Wolf-Rayet stars and 
supergiant B[e] stars. 
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Fig. 5: OGLE photometry. Phases are computed with re¬ 
spect to the photometric circular-orbit Tq values and peri¬ 
ods reported in Table OGLE-HI and OGLE-IV mag¬ 
nitudes are shown in grey and black, respectively; the 
OGLE-HI results for VFTS 450 have been offset by —CP 174 
(cf. |2.2| ). The V- and Ic-band measurements were not 
quasi-simultaneous, and OGLE-HI colours have been com¬ 
puted from phase-binned Ic results, which are ~10x as 
numerous as the V measurements. 


3. Spectroscopic orbits 

3.1. Radial-velocity measurements 

3.1.1. Primaries 

Radial velocities for the primary components were rea¬ 
sonably straightforward to measure using relatively un¬ 
blended Sim and Hei lines. We used the results of the 
model-atmosphere analyses reported in Section [^ to iden¬ 
tify suitable tlusty models to employ as templates in a 
cross-correlation analysis. Results, which are insensitive to 
the precise choice of model template, are incorporated into 
Tables and 1^ (on-line); the dispersion in velocities from 
different lines, and residuals from the orbital solutions dis¬ 
cussed in Section [3.2[ are consistent with measurement er¬ 
rors of <10 km s“-^. 


4 


colour 























I.D. Howarth et al.: Two massive binaries 




Wavelength (A) 


Fig. 6: Selected helium lines in the spectrum of VFTS 450 
near quadrature; smooth red curves show gaussian fits to 
the data. Although these helium lines give consistent results 
for the primary (at ^+400 km s~ ^), the secondary veloci¬ 
ties are discordant (Section 3.1.2 for reference, the dashed 
vertical line indicates the measured secondary He i A4471 
velocity). This figure also illustrates the differen ces in 
He i:He ii line ratios in the two components (Section |4.2[). 


3.1.2. Secondaries 

In order to measure the much weaker secondary spectra, 
we merged LR02/LR03 spectra taken on any given night 
using a weighted mean, with a sigma-clipping algorithm to 
exclude cosmic-ray events and other flaws. (Multiple obser¬ 
vations taken at a given spectrograph setting on any one 
night span <1% of the orbital periods that we report in 
Table ^ and may therefore be combined without special 
procedures to compensate for binary motion.) Uncertain 
corrections for echelle blaze render measurements in the 
UVES spectra unreliable. 

The secondary velocities were measured by direct fitting 
of gaussians, but the shallowness and breadth of the lines 
make the results quite sensitive to the adopted rectification. 
Repeat measurements and residuals to model fits are both 
consistent with typical measurement errors of ~35 km s“^. 


VFTS 450 Because of blending with features in the pri¬ 
mary spectrum, we did not attempt radial-velocity mea¬ 
surements of the secondary at phases near conjunctions 
(100 < Vprimary ^ 400 km s“^). The helium lines in the 
secondary spectrum show poor agreement, as illustrated in 
Fig.|6) the He ii A4200 velocities are generally - though not 
consistently - some ~100 km s“^ more positive than found 
for Hell A4541 or Hei A4471. Given the shallowness of the 
lines in the secondary spectra, we cannot rule out that rec¬ 
tification difficulties contribute to this problem. In practice, 
we rely principally on results for A4200, which gives consis¬ 
tent results and which is not subject to significant blending 
(cp., e.g., secondary A4541, which can be affected by Sim 
A4552 in the primary spectrum). 


VFTS 652 The absence of He ii lines in the primary spec¬ 
trum renders measurement of the weak He ii A4541 line in 
the secondary reasonably straightforward in both LR02 and 


Fig. 7: VFT S 652 as a double-lined spectroscopic binary 
(Section|4^. Wavelengths of selected lines are shown in the 
rest frame of the system centre of mass, together with the 
observed orbital displacements of N ii A4530, Si ill A4552 
(primary spectrum), and He ii A4541 (secondary spectrum). 



Fig. 8: Spectroscopic orbit for VFTS 450; orbital phases re¬ 
fer to the circular-orbit Tg (from Table[^, as do the (O—C) 
residuals for the primary shown at the oottom of the plot. 
The eccentric-orbit solution for the primary is shown as a 
dotted line (which may appear as a continuous grey line if 
viewed at low resolution), and the circular orbit for the sec¬ 
ondary as a dash-dot line. Diamonds show the He ii A4200 
velocities measured in the secondary’s spectrum, and used 
to estimate the mass ratio. 


LR03 spectra. He ii A4200 gives consistent, but somewhat 
more scattered, results. 


3.2. Results 

The primaries’ spectroscopic orbits are summarized in 
Table and are illustrated in Figs. and We adopted 
uniform weighting for our final solutions, but other weight¬ 
ing schemes result in unimportant changes to the orbital 
parameter s. According to the form ulation of the F test de¬ 
scribed by Lucy & Sweeney (1971), the orbital eccentricities 
are formally significant with > 99% confidence. However, 
the apparent eccentricities are quite small, and we caution 
that they may not reflect the true centre-of-mass motions. 
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Table 5: Radial-velocity orbital solutions, for circular and eccentric orbits. 


Parameter 

VFTS 450 
circ. ecc. 

VFTS 652 
circ. ecc. 

Pphot (d) 

6.892583 

8.589555 

rb 

0.000039 

0.000090 

TqP (MJD) 

54761.421 


54665.242 


± 

0.007 


0.014 


^spec (d) 

6.892325 

6.892242 

8.589534 

8.589413 

zb 

0.000124 

0.000080 

0.000089 

0.000067 

7 (km s“^) 

248.46 

246.26 

253.09 

255.32 

zb 

1.61 

1.10 

0.83 

0.65 

K (km s"^) 

208.8 

207.6 

201.6 

199.8 

± 

2.2 

1.5 

1.4 

1.0 

e 

= 0 

0.0722 

=0 

0.0443 

± 


0.0071 


0.0050 



356.2 


10.0 

zb 


6.3 


5.4 

To (MJD) 

54761.268 

54761.214 

54897.050 

54897.300 

zb 

0.021 

0.116 

0.011 

0.127 

f{M) (Mo) 

6.51 

6.35 

7.31 

7.10 

zb 

0.20 

0.13 

0.15 

0.11 

ai sini (Ro) 

28.43 

28.20 

34.22 

33.88 

zb 

0.29 

0.20 

0.24 

0.18 

r.m.s. residual (km s“^) 

13.8 

9.1 

7.5 

5.4 

q{=MilM2) 


0.61 


0.40 

zb 


0.05 


0.05 

Misin^(i) (Mq) 

10.4 

10.2 

5.6 

5.5 

zb 

1.5 

1.4 

1.2 

1.2 

M'2sin^(i) (Mq) 

16.9 

16.5 

14.2 

13.8 

± 

1.2 

1.0 

1.1 

1.1 


Notes. Parameters are based on primary-star radial velocities, excepting the photometric period, Pphot; the photometrically 
determined value of the time of circular-orbit maximum radial velocity, Tq; and the mass ratio q (Section 3.2 I. Note that the 
To parameter has different meanings for circular and eccentric orbits (times of maximum velocity and of periastron passage, 
respectively); for both targets, the numerical values from the spectroscopic solutions are coincidentally similar only because ~ 0° 
in each case. 



Orbital Phase 

Fig. 9: Spectroscopic orbit for VFTS 652; details as for 
Fig.lH except that He ii A 4541 velocities are shown for 
the secondary. 



Given the considerable uncertainties in the radial- 
velocity measurements of the secondaries, we chose a sim¬ 
ple but robust method to estimate the mass ratio for each 
system, namely, a linear regressio n o f the secondary veloc¬ 
ities on the primary values (Fig. 101. The gradient yields 


Fig. 10: Spectroscopic mass-ratio determination; secondary 
velocities have been offset as indicated for display purposes 
(primary velocities unchanged). The slope of the linear fit 
gives the mass ratio directly. 
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Fig. 11: Disentangled component spectra; the y-axis scaling is approximate (and depends on the primary:secondary 
continuum flux ratios). The Balmer lines and, to a lesser extent, He i are corrupted by nebular emission (marked as ‘n’); 
broad features, such as the A4430 diffuse interstellar band, are rectified out. 


the mass ratio directly, independently of all other param¬ 
eters; results are included in Table The means of the 
observed minus predicted secondary orbital velocities (i.e., 
the differences between primary and secondary 7 velocities) 
are -bl0.0±8.5 (s.e.) and -b29.5±6.5 km s"! for VFTS 450 
and 652, respectively; differences in 7 velocities such as that 
shown by VFTS 652 h ave been found previously in ‘windy’ 
massive binar ies (e.g., Niemela fc Morr^|1986[ Niemela & 
Bassino|[l994 |. 


Although the spectroscopic period determined for 
VFTS 450 differs from the photometric value by ^1.7-cr, 
we don’t consider this to be significant evidence for period 
changes, given that the OGLE-III and OGLE-IV datasets 
are in good mutual agreement, and span the spectroscopic 
epochs. 


3.3. Disentangling 


In principle, an alternative approach to the spectroscopic- 
orbit modelling is a simultaneous solution of the individual 
component spectra and the orbital characteristics (‘disen¬ 
tangling’; cf., e.g., Hadrava|2004 ). Because of the weakness 
of the secondary spectra we instead chose the simpler op¬ 
tion of reconstructing the separate component spectra in 
the more extensive LR02 datasets from the ‘known’ veloci¬ 


ties, using CRES (|Ilijic|2004 ). We explored the consequences 
of using observed velocities or those calculated from the or- 


: using 

bital solution, including using mass ratios in the range 0.5- 
1.0 when computing the secondary spectra. We found the 


results to be quite robust to these factors (the corollary be¬ 
ing that the technique cannot recover a precise mass ratio 
for these data). The resulting spectra of individual com¬ 
ponents are shown in Fig. [m they have better S:N than 
any individual spectrum, but the y scaling is arbitrary, and 
nebular emission contaminates some key lines. 


4. Model-atmosphere analysis 

4.1. Methodology 

Model-atmosphere analyses of both systems were per¬ 
formed using a grid of line-blanketed non-LTE tlusty 


models at LMG metallicity ( 

Hubeny|1988 

Hubeny & Lanz 

19951 IHubeny et al. 

19981 for more details of the grid see 

Ryans et al. 2003 

Dufton et al. 

x? -:—r: 

2005 

). The analyses as- 


acterized by a single set of atmospheric parameters, and 
that hydrostatic, plane-parallel structures are appropriate. 
Depending on the adequacy or otherwise of these assump¬ 
tions, the results may be subject to significant (and largely 
unquantifiable) systematic errors, and should therefore be 
interpreted with due caution. 

For the primary spectra, the atmospheric parameters 
were estimated from the Sim and Siiv line strengths, to¬ 
gether with the H i and He ii profiles. The lower quality of 
the secondary spectra allowed only relatively rough esti¬ 
mates of parameter values to be made, using the Hi and 
He II lines; the microturbulence was indeterminate (and 
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Table 6 : Equivalent widths of the primary in the integrated 


spectra; cf. 14.1.1 


Species 

Wavelength 

(A) 

Region 

IUa ( 
[450] 

mA) 

[652] 

Nil 

3995.0 

LR02 

26 

185 

Nil 

4227.7 

LR02 


39 

Nil 

4447.0 

LR02 


70 

Nil 

4601.5 

LR03 


113 

Nil 

4607.2 

LR03 


117 

Nil 

4613.9 

LR03 


77 

Nil 

4621.4 

LR03 


84 

Nil 

4630.5 

LR03 

41 

189 

Nil 

4774.2 

LR03 


12 

Nil 

4788.1 

LR03 


24 

Nil 

4803.3 

LR03 


51 

Nil 

4994.4 

LR03 


40 

On 

4185.4 

LR02 


20 

On 

4317.0 

LR02 


61 

On 

4319.6 

LR02 


59 

On 

4366.9 

LR02 


63 

On 

4393.9 

LR02 


14 

On 

4414.9 

LR02 

26 

102 

On 

4417.0 

LR02 

16 

72 

On 

4452.4 

LR02 


12 

On 

4591.0 

LR03 

24 

62 

On 

4596.0 

LR03 


44 

On 

4661.6 

LR03 

34 

96 

Mg II 

4481.0 

LR02 

40 

70 

Sim 

4552.6 

LR02 

98 

230 

Sim 


LR03 

84 

245 

Sim 

4567.8 

LR03 

69 

198 

Sim 

4574.8 

LR03 

39 

116 

Siiv 

4088.9 

LR02 

322 

104 

Si IV 

4116.1 

LR02 

233 

84 


unimportant) for these lines, and we assumed appropriate 
values. 

The analyses were based principally on the disentan¬ 
gled spectra, which have better signahnoise ratios than any 
individual spectrum (and, of course, should be free from 
blending), although cross-checks were made against the di¬ 
rectly observed spectra, particularly for the Balmer lines, 
which suffer nebular contamination. 

A complication is the uncertainty in the relative flux 
contributions of the two components in each system; with¬ 
out ancillary information, it is impossible to distinguish 
between a strong continuum with weak lines and a weak 
continuum with strong lines. We addressed this issue 
by supposing that the primary [secondary] contributes a 
wavelength-independent fraction [J^ 2 ,= 1 — -Ai] of the 
rectified continuum flux, and adjusted this fraction as nec¬ 
essary. 


4.1.1. Metal-line equivalent widths 

The data quality allows abundance analyses to be con¬ 
ducted for the primary spectra. For this purpose, equivalent 
widths were measured by fitting theoretical profiles to the 
observations, using a least-squares technique. In the LR03 
region directly observed spectra were used; for the LR02 
setting, the disentangled spectra were employed (to take 
advantage of the improved S/N), and the results scaled to 


recover the W\ values that would be determined in the di¬ 
rectly observed spectra. Results are summarized in Table 

The A4552 Sim line falls in the region of overlap be¬ 
tween LR02 and LR03 spectra; results from both spectro¬ 
graph configurations are separately listed in the Table. The 
values from the two settings differ by ~15mA, in opposite 
senses for VETS 450 and VETS 652. This is probably a fair 
reflection of observational uncertainties (including, poten¬ 
tially, temporal/orbital variations, although comparison of 
LR02 spectra taken at different epochs shows no evidence 
for substantial changes in line strengths). 

The disentangling results show that the absorption lines 
of metals can be safely attributed to the primary spectra 
(even in the LR03 data), but of course the measurements in 
Table 1^ have to be scaled by the appropriate value when 
performing an abundance analysis. 

4.2. VFTS450 

The LR02 primary spectrum shows hydrogen and neutral & 
ionized helium lines, together with strong metal lines (par¬ 
ticularly On and Siiv). By contrast, only the hydrogen 
and ionized helium lines are clearly seen in the secondary 
spectrum, with the former being badly contaminated by 
nebular emission. The Hei lines at 4387, 4471A are prob¬ 
ably present, with less convincing evidence for Siiv A4089 
and Niii AA4097, 4510-4534. 


Primary: The primary spectrum was first analysed by as¬ 
suming no secondary contamination (i.e., J -2 = 0 ), leading 
to estimates of and logg that are, in practice, lower 
limits to allowable values. To investigate the sensitivity of 
the results to spectral contamination by the secondary com¬ 
ponent, the analysis was repeated for J '2 = 0.25, 0.5, span¬ 
ning the range of plausible values; results are summarized 
in Table [T] 

Gravities were estimated from the Hi5 and H 7 line pro¬ 
files, with results agreeing to better than 0.1 dex. The ef¬ 
fective temperature was estimated from the He ii spectrum 
(by assuming a normal helium abundance), as results from 
the silicon ionization equilibrium were found to be sensitive 
to the microturbulence, and were used to determine that 
parameter. (Because of the relatively large value of the 


Sim triplet lines at 4552-4574A lie near the linear part of 
the curve of growth, and hence are not particularly sensitive 
to the microturbulence.) 

Given ^ and iF, element abundances were estimated 
from the equivalent widths listed in Table with stan¬ 
dard deviations of 0 . 1 - 0.2 dex implied by the individual 
oxygen estimates. To obtain approximately ‘normal’ LMG 
abundances for magnesium and silicon requires F 2 — 0.25, 
which represents our ‘best-bet’ model. There then appears 
to be a significant surface-nitrogen enhancement approach¬ 
ing 1 dex, and an oxygen depletion of ^0.3 dex. 

The range of parameter estimates, together with the 
agreement between theoretical and observed profiles, leads 
us to adopt modelling uncertainties of ±1 kK in 
±0.2 dex in logg, and 2 km s“^ in Uncertainties on the 
abundances are difficult to address precisely; the adopted 
uncertainties in the atmospheric parameters alone translate 
into typical errors of 0.1 5 dex for both nitrogen and oxygen, 
(see Hunter et al.||2007 for more details). Varying the sec- 
ondary contribution (± 2 ) contributes significant additional 
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Table 7: Summary of stellar-atmosphere parameter determinations. 


Star 

±■1 

T 2 

Teff 

logs 


Vq sin i 


Abundances 





(kK) 

(cgs) 

(km s ’^) 

(kms-Q 

N 

0 

Mg 

Si 

N/0 

VFTS 450 p 

1.00 

0.00 

25.5 

2.7 

15 

99 

7.32 

7.75 

6.93 

6.70 

-0.43 

*VFTS 450 p 

0.75 

0.25 

27.0 

2.9 

15 

99 

7.63 

8.06 

7.13 

7.17 

-0.43 

VFTS 450 p 

0.50 

0.50 

29.0 

3.1 

18 

99 

8.02 

8.48 

7.34 

7.68 

-0.46 

*VFTS 450 s 

0.75 

0.25 

33.5:35.0 

<3.8 

[10] 

320 






VFTS 450 s 

0.50 

0.50 

28.5:31.0 

<3.4 

[10] 

320 






VFTS 450 s 

0.00 

1.00 

25.5:27.5 

<3.0 

[10] 

320 






VFTS 652 p 

1.00 

0.00 

21.3 

2.6 

9 

83 

8.04 

7.88 

6.89 

6.94 

±0.16 

*VFTS 652 p 

0.75 

0.25 

22.2 

2.8 

13 

83 

8.13 

7.96 

7.00 

7.07 

±0.17 

VFTS 652 p 

0.50 

0.50 

23.0 

3.0 

19 

83 

8.31 

8.08 

7.15 

7.28 

±0.23 

*VFTS 652 s 

0.75 

0.25 

35.0 

3.7 

[5] 

260 






VFTS 652 s 

0.50 

0.50 

30.0 

3.1 

[10] 

260 






VFTS 652 s 

0.00 

1.00 

26.0 

2.7 

[10] 

260 









Reference LMC baseline abundance: 

6.90 

8.35 

7.05 

7.20 

-1.45 


Notes. T\.i are the adopted fractional contributions of the primary, secondary (p, s) components to the total continnum flux, and 


^ is the microtnrbulence (with assumed values given in square brackets). Baseline LMC abundances are taken from Hunter et al. 
(20071 on a logarithmic scale where the abundance of hydrogen by number = 12). Abundances are given to two decimal places to 
clarify the sensitivity t o T values, and not to indicate the accuracy of the determinations (for which reali stic errors are > 0.2 dex, 
as discussed in Section [4.4[ ). Similarly, projected rotation velocities are probably good to only ~10-20% ( ^4.5| . Preferred solutions 
are indicated by asterisks. 


uncertainty to the absolute abundances, but has little effect 
on the N:0 abundance ratio; the inference of a significant 
surface-nitrogen enhancement appears to be robust. 

Secondary: The weakness of the secondary’s absorption 
lines makes an atmospheric analysis difficult, and our re¬ 
sults should be treated with caution. Nonetheless, the 
Hell absorption lines do provide useful diagnostics, princi¬ 
pally for the effective temperature. Additionally, although 
the Balmer-series lines are badly contaminated by nebn- 
lar emission, the lack of significant Stark-broadened wings 
sets an upper limit on the surface gravity. As an exercise 
in defining the range of possible parameter space, we con¬ 
ducted analyses assuming that the continuum was entirely 
due to the secondary {J ^2 = 1 ), along with two cases con¬ 
sidered for the primary {T 2 = 0.25 and 0.5); the results 
are summarized in Table Note that the gravity limits are 
appropriate for the upper limit of the effective-temperature 
range - a lower effective temperature would lead to lower 
gravity limit. 

4.3. VFTS 652 

The spectra of VFTS 652 show a rich metal-line spectrum 
for the primary, in accord with its classification as a B- 
type supergiant. The secondary spectrum shows convincing 
evidence for the presence of Hell lines (Fig. [^, together 
with Stark-broadened wings in the Balmer series; He i lines 
also appear to be present. 

Primary: As for VFTS 450, we evaluated parameters for 
= 0, 0.25, and 0.5. The effective temperature was esti¬ 
mated from the silicon ionization equilibrium, the gravity 
from the Balmer-line profiles, and the microturbulence from 
the relative strengths of lines in the Si ill triplet. Results are 
summarized in Table 0 


To obtain ‘normal’ LMC abundances for magnesium 
and silicon requires J -2 = 0.25-0.50. Abundances were also 
derived for nitrogen and oxygen, with the scatter among es¬ 
timates from individual lines being in the range 0 . 1 - 0.2 dex. 
Using the same criteria as for VFTS 450 leads to estimated 
uncertainties of ±1 kK, ±0.1 dex, and ±3 km s“^, for Tes, 
log 5 , and ^ respectively. These in turn imply uncertain¬ 
ties of typically 0 . 2 - 0 .3 dex for the nitrogen and oxygen 
abundances (but significantly less for the N:0 abundance 
ratio). The effects of varying the dilution factor are rela¬ 
tively small, leading to adopted final errors of 0.3 dex for 
these elements; a conservative error estimate on the N:0 
ratio is ^ 0.2 dex. 

At J ^2 = 0.25, the oxygen abundance is possibly under- 
abundant, by ^0.4 dex compared to the LMC baseline, 
while nitrogen is again clearly enhanced, by more than 
1.0 dex. Qualitatively similar conclusions follow for other 
dilution factors; regardless of the dilution factor adopted, 
the N:0 abundance ratio is 1.6-1.7 dex higher than for the 
adopted LMC baseline abundances. 


Secondary: Again as an exercise in defining the range of 
possible parameter space, we used the Hen AA4200, 4541 
and available Balmer lines to estimate atmospheric param¬ 
eters for several dilution factors. The results are summa¬ 
rized in Table 0 for T 2 = 0.25, inferred stellar parameters 
lie at the boundary of, or just outside, our grid of tlusty 
models. Representative values of the microturbulence were 
adopted but varying these by reasonable amounts would 
have a negligible effect on our estimates. 

The absence of detectable metal lines in the secondary 
spectrum is consistent with the adopted parameters, sec¬ 
ondary flux level, and signaknoise ratio. 
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4.4. Uncertainties 

Atmospheric parameters: The primary spectrum of 
VFTS 652 was easier to model than that of VFTS 450, 
with better internal agreement between different spectral 
features. Flowever, for both systems the inclusion of 
dilution by the secondary leads to only relatively small 
changes in fit parameters. We therefore expect the error 
estimates for the primaries’ and logg values discussed 
above to be reasonable. 

Parameter estimates for the secondary components are 
considerably less secure. The Teff and log g values estimated 
from the Hen and Hi profiles are moderately sensitive to 
the choice of the dilution factor, 7 ' 2 - It is therefore dif¬ 
ficult to assess realistic error estimates, although in gen¬ 
eral terms we consider the T^g values to be more reli¬ 
able than those for log 5 (reflecting the greater reliability 
of the deconvolved Hell spectrum); for the ‘worst case’ of 
VFTS 450, we estimate an uncertainty in Teff of perhaps 
±4 kK. Nonetheless, provided that the secondary spectra 
are formed in the photospheres of the secondary stars, it 
seems secure that the primary is the cooler component in 
each system. 


Dilution factor: From the general characteristics of the spec¬ 
tra, we are confident that 7^2 is certainly less than 0.5 
for each system. Unfortunately, the magnesium and sili¬ 
con abundances used to constrain on the dilution factor 
also depend the atmospheric parameters (and hence do not 
provide particularly strong limits). However, that the sec¬ 
ondary spectra can be measured at all implies 7-2 > 0 . 1 . 
Hence dilution factors of 7^2 — 0.25 ± 0.1 would appear to 
be reasonable for both systems. 


Abundance estimates: These are available only for the pri¬ 
mary components, and are quoted in Tablej^to two decimal 
places in order to illustrate the sensitivity to the adopted 
dilution factors (and not to indicate their accuracy). The 
most striking results are the enhanced surface-nitrogen 
abundances (0.7dex for VFTS 450, 1.2 dex for VFTS 652) 
and nitrogen:oxygen abundance ratios ( 1.0 dex and 1.6 re¬ 
spectively); although there are additional uncertainties as¬ 
sociated with binarity, the relatively normal abundance es¬ 
timates for other elements strongly support large nitrogen 
enhancements in both stars. 


The situation is less clear for oxygen, with implied un¬ 
derabundances of 0.3-0.5 dex, compared to uncertainties of 
± 0 . 2 - 0 .3 dex; such underabundances are, however, consis¬ 
tent with those predicted from LMC single-star evolution¬ 
ary mode ls that yield nitrog en enhancements of a factor 
^10 (e.g., Brott et ar]|2011a ). 


spectra ((3.31, but checks were performed using the directly 
observed data. 

For the primaries, we used N ii AA3995, 4447, 
Sim A4552, Siiv A4089, and Siii A4253. Adopted UgSini 
values are averages of results from all lines in each target; 
although the dispersion from different lines is <5%, system¬ 


atic effects can be important (see notes in Sundqvist et al. 
2013 Simon-Di'az fc Herrero|2014 1, and realistic uncertain¬ 


ties are perhaps ^ 10 - 20 %, following arguments given by 
(20151. For the secon^ries we were limited 
e He II lines (Section]^, but systematic ef- 


McEvoy et al. 
to the availab’ 


fects should be negligible, and likely uncertainties are on 
the order of ~ ± 10 %. 

Results are included in Tablej^ The primaries’ rotation 
velocities are rather high when considered in the context 
of the VFTS sample of single late-O/early-B supergiants 
(McEvoy et al. 2015), but a re close to values expected for 
synchronous rotation ( |5.2[ ); however, the secondaries’ ro¬ 
tations are exceptionally rapid. 

Estimates of v sin i from simple profile fitting of the pri¬ 
mary spectra are 105 ± 6 and 91 ± 4 km s“^ (VETS 450, 
652), ~5-10% larger than VeSini measurements from the 
Eourier Transform methodology; however, the estimates are 
consistent within the uncertainties. Thus while there may 
be a macroturbulent contribution to line broadening in the 
primary spectra, its extent is difficult to quantify usefully. 
Rotational broadening dominates the secondary spectra, 
rendering estimates of any macroturbulence contribution 
impossible. 


5. Discussion 


The secondary-spectrum radial velocities indicate that the 
spectroscopically less conspicuous star is, in each system, 
the more massive component. This conclusion is supported 


by photometric considerations (§(5.1.1 

5.1.2 

), and so ap- 

pears to be a robust conclusion, even if 1 

he secondary spec- 


trum is only an approximate tracer of the secondary star’s 
centre-of-mass motions. 

Taken at face value, the spectroscopic analysis also in¬ 
dicates the secondaries to be the hotter components and 
hence, being fainter, the smaller. However, in each system, 
the secondary spectrum has the larger UgSinf; we cannot, 
therefore, assume corotation in order to constrain radii or 
inclinations. In principle, a light-curve analysis can yield 
this information, but eclipses, if they occur at all, are very 
shallow (Eig. [^, leading to poorly constrained solutions, 
and our initial attempts in this direction have yielded un¬ 
physical results. The lack of deep eclipses, coupled with sig¬ 
nificant ‘ellipsoidal’ variations, does, though, indicate both 
that the primaries are close to filling their Roche lobes, 
and that the systems are observed at intermediate orbital 
inclinations. 


4.5. Projected rotation velocities 

Estimates of the projected equatorial rotation veloci¬ 
ties, UgSini, were obtained by using a Fourier tech¬ 
nique similar to that adopted in othe r VFTS rotational- 
velocity investigations (see, fo r example Dufton et al.|201^ 


Ramfrez-Agudelo et al. 12013), and by simple profile fitting, cf., e.g. 


5.1. System constraints 

The systems’ absolute magnitudes can be determined from 
the apparent magnitudes (we use the mean V,Ic OGLE 
results f rom Table [4|), t he LMC distance modu lus (18.5; 
Schaeferl 120081 


which yields the line-width parameter v sin i (which has 
contributions from both rotation and macroturbulence). 
Measurements were principally made on the disentangled 


Pietrzyiiski et al. 2013), and the 


reddenings. 

We estimate {B — V)q and {V — Ic)o by taking flux- 
weighted averages of empirical intrinsic colours for each 
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comp onent as a function of spectral type (from Wegner 
19941, and of model colours as a function of tempera¬ 


ture (making use of synthetic photometry from the LMC- 


abundance Atlas models reported by Howarth 


2011 ). 


The two sources of intrinsic colours, wit h two observed 
colours (Table 0, and a reddening law (Howarth 1983 


A{Ic)/E{B — V) = 1.84), then yield four separate esti¬ 
mates of E{B — V) and E{V — Ic), whence four estimates of 
M{y) and M{Ic) for each system. The dispersions in these 
estimates are small, and we simply adopt mean values; for 
quantitative results we rely principally on the /c-band re¬ 
sults, since the extinction and the sensitivity of flux to tem¬ 
perature are both slightly less here than at V. (Adopting 
the absolute V magnitudes introduces only minor changes 


to the numerical results, as illustrated in Fig. 13 


Absolute magnitudes for the individual components fol¬ 
low from the continuum flux ratio, E 2 /J-i[= ~ ■^ 2 )]- 

Coupling these with the corresponding surface fluxes (from 
model atmospheres at the spectroscopically-determined ef¬ 
fective temperatures) gives the stellar radii. 

The observed Oi sin i value, together with the mass ratio 
g, gives both the projected semi-major axis asini, and the 
projected Roche-lobe radii, i?L(l, 2) sini (convenie ntly eval¬ 
uated using the analytical approximation given by|Eggleton| 


1983). Requiring the primary’s radius not to exceed its 
Roche-lobe radius sets a limit on sin i for a given q\ or, alter¬ 
natively, limits possible values for q (by setting sini = 1 ). 

With values for i?, q, and i in hand, other para¬ 
meters (L, M, etc.) follow straightforwardly, given the pri¬ 
mary’s spectroscopic orbit. 


{B - V) baseline [E{B - V) = 0.40 ^ E{V - Ic) = 0.50] 
than for {V - Ic) [E{B -V) = 0A7 ^ EjV - Ic) = 0.59]. 

The same reasoning as applied in §5.1.1| gives an up¬ 
per limit on the primary’s radius of i?i < 25.SR© (±0.9, 
±1.2Rq). This implies q < 0.98±0.05 (again indicating that 
the secondary is the more massive component), or, adopting 
the spectroscopic mass ratio, sinz < 0.66; while the spectro¬ 
scopic brightness ratio of ~3:1 implies i?i ^2 — 22.1,8.5 R©. 


5.2. A first estimate of system parameters from spectroscopy 

The substantial photometric variability strongly suggests 
that the primaries fill, or very nearly fill, their Roche lobes, 
as do the various indicators of lobe-overflow mass trans¬ 


fer, discussed further below (]5.5). With this assu mption, 
and using the procedures outlined in Section |5.11 we can 
make a first estimate of approximate actual system param¬ 
eters, which are summarized in Table (columns headed 
‘Ml’). The Table also explores the sensitivity of derived 
quantities to input parameters. Masses are the least well de¬ 
termined variables, principally because of the third-power 
dependence on sinz. 

The inferred inclinations are consistent with the ab¬ 
sence of clear eclipses, and there is tolerable agreement be¬ 
tween the orbital and model-atmosphere estimates of log 5 , 
although the spectroscopic determinations are ~0.3 dex 
smaller]^ For these first-pass parameter estimates, pro¬ 
jected equatorial corotation velocities are in good agree¬ 
ment with the primaries’ observed values, but the secon¬ 
daries appear to be rotating considerably faster than syn¬ 
chronous (although well below critical). 


5.1.1. VFTS 450 

We find M(y) = -6.35 ± 0.11, M(/c) = -6.01 ± 0.07, 
where the errors are standard deviations of the four in¬ 
dividual estimates (which are not independent). The in¬ 
ferred reddening is slightly larger for the {B — V) baseline 
[E{B -V) = 0.45 ^ E{V - Ic) = 0.57] than it is for 
[V - Ic) [EiB -V)= 0.39 ^ E(y - Ic) = 0.49]. 

The upper limit on the primary radius, assuming that it 
contributes all the /c-band light, is 25.4 R©, for Teff(l) = 
27 kK (±0.8, ±1.1R© for ATes = ±1 kK, AM = TO.p^. 
More realistically, using the spectroscopic (^ H-bancQ 
brightness ratio of ~3:1, the implied radii are i?i _2 — 22!o, 
10.1 R©, with uncertainties on the order of 10%. 

For a primary radius Ri < 25.4±1.0 R© we find 
q < 1.24 ±0.06. This can be considered a rather firm upper 
limit, as it depends only on the absolute magnitude, the pri¬ 
mary’s effective temperature, and its radial-velocity curve, 
all of which are reasonably well established; this analysis 
therefore suggests that the secondary is very probably the 
more massive component (independently of the secondary 
radial-velocity curve). Adopting the spectroscopic mass ra¬ 
tio of 0.61 implies sinz < 0.70, where the equality corre¬ 
sponds to a lobe-filling primary. 


5.1.2. VFTS 652 

We estimate M{V) = -5.86±0.14, M(/c) = -5.67±0.08; 
in this case, the inferred reddening is slightly smaller for the 

^ For completeness, we adjust the ±2 values as a function of 
wavelength by using model-atmosphere fluxes, though this has 
negligible consequences. 


5.3. Light-curves 

Model light-curves for parameters in the region of the ‘Ml’ 
solutions fail to reproduce the amplitudes of the observed 
light-curves. The observed ‘ellipsoidal’ variations (Fig. 
imply that the primary in each system must be very close 
to filling its Robe lobe, but if the secondary is hotter and 
fainter than the primary, while being more massive, then 
it must significantly underfill its Roche lobe, regardless 
of detailed numerical parameter values. The amplitude of 
orbital photometric variability under these circumstances 
does not substantially exceed ^ 0±1 over a range of mass 
ratios and inclinations - about half the observed ampli¬ 
tudes. Varying the spectroscopically inferred parameters 
over plausible ranges cannot overcome this discrepancy; the 
only way to reproduce the light-curve amplitude by conven¬ 
tional models is to adopt an overcontact (or double-contact) 
configuration, but this would imply spectroscopically more 
conspicuous secondaries. 


5.4. Evolutionary considerations 


The schematic ‘Ml’ system parameters are plotted in an 
H-R diagram in Fig. |12[ evolutionary tracks for non¬ 
rotating single stars at LMC metallicity, from|Brott et al.| 
( |2011a|bD , are also shown. This figure discloses a further 
problem: although the dynamical masses estimated for the 
primaries (i.e., the cooler, less massive, lobe-filling compo¬ 
nents) are in reasonably good agreement with the single- 


^ Corrections for centrifugal forces, Alogg ~ (ue sinz)^/i?,, 
are <0.03. 
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Table 8: System parameters based on spectroscopy. 
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Log 


Fig. 12: Stellar parameters from Table plotted in the 
Hertzprung-Russell diagram. Error bars illustrate uncer¬ 
tainties of ±1 kK on primaries (larger symbols) and ±2 kK 
on secondaries, and the sums in quadrature of the error 
ranges on L listed in Table Dynamical masses are indi¬ 
cated in squa re brackets. Evolutionary tracks from |Brott| 


et al. (2011a I for single, non-rotating stars are shown for 
comparison, labelled by ZAMS mass. 


star tracks, the secondaries are significantly under-luminous 
for their dynamical masses - and standard binary evolution 
cannot produce this outcome. 

Nevertheless, it is unlikely that these systems can be 
anything other than hot, massive counterparts of typical 
Algol-type binaries, in the slow (nuclear-timescale) phase 
of Case A mass transfer. Mass transfer (including common- 
envelope evolution) in a more evolved configuration would 
produce a helium star (a Wolf-Rayet star at these masses), 
which would be spectroscopically conspicuous. Chemically 
homogeneous evolution of the rapidly rotating secondaries 
would lead to significantly higher effective temperatures, 
and can probably also be excluded. 


5.5. Parameter-space exploration 


Given the difficulties encountered in reconciling spectro¬ 
scopic, photometric, and evolutionary constraints, the fin¬ 
ger of doubt points most directly at the unqualified attri¬ 
bution of the observed secondary spectra to the secondary 
stars’ photospher es. T aken together, the breadth of the ab¬ 
sorption profiles (j|4.5 ), the discrepancies in radial velocities 
from different lines ((3.1.2), the anomalous dou ble-p eaked 
Ha profiles (Fig. [^, unusual near-IR colours ((2.2), and 
the general elusiveness of the secondary spectra, all suggest 
the possibility that the secondary’s photospheric spectrum 
in each system may be modified, or even concealed, by an 
accretion disk (which would also be consistent with lobe¬ 
filling primaries). 

These observ ed properties are reminiscen t of the W Ser 
class of binaries ( Plavec|1980 Tarasov|2000 1 . Although the 
VFTS targets have higher masses and shorter periods than 
is typical for this group, their general characteristics, in¬ 
cluding lobe-filling, synchronously rotating primaries, and 


IR excesses, are in accord with this notion ([Andersen fc 


Nordstrom |1989 Mennickent & Kolaczkowski PUlOp , and 
there are clear amilarities to related early-type systems 


such as RY Set an d V453 Sco (|Grundstrom et al. 2007 
Josephs et al.||200T ). 


The VFTS binaries studied here may well, therefore, 
have secondary spectra that are contaminated by, or arise 
in, accretion disks. In recognition of this possibility, we ex¬ 
plore a broader parameter space for the systems, and par¬ 
ticularly for the secondary components. 


5.5.1. Mass ratios, brightness ratios 


For heuristic purposes, we first consider the consequences 
of adopting mass ratios q and continuum brightness ratios 
in ranges outside those directly inferred from the 
secondary spectra. As described in Section |5.H a bright¬ 
ness ratio and component temperatures yield the stellar 
radii; given the primary radius, a mass ratio gives the in¬ 
clination (for a lobe-filling primary), and hence, from the 
spectroscopic orbit, the masses. Luminosities follow from 
the radii and temperatures. 

The basic inputs we adopt for each system are therefore 
(i) the absolute magnitude; (ii) the orbital period; (iii) the 
primary’s orbital-velocity semi-amplitude; and (iv) the two 
components’ effective temperatures. With fixed values for 
these quantities, the system characteristics are fully spec¬ 
ified by q and T^jTx (assuming that the primary fills its 
Roche lobe). 

Results in the mass-luminosity plane are illustrated in 
Fig. In this figure, M-L curves are shown for selected 
specific values of the mass ratio g, over a range in 
The curves are insensitive to q for the primary stars, but 
not for the secondaries. 

Also shown in the Figure are lines of constant orbital 
inclination. For a given q, then a particular inclination cor¬ 
responds to a specific this {g, TijJ-x) pair yields 

the full set of other parameters, including M and L. Thus 
a constant inclination corresponds to a curve in the M- 
L plane (again, assuming that the primary fills its Roche 
lobe). 

Other than at advanced evolutionary stages, a single 
hot star will norr ually lie between the ZAMS and TAMS 
M-L loci in Fig. 131"^ If binary evolution is to produce a 
secondary that is not underluminous for its mass, then this 
component must lie somewhere above the ZAMS M-L lo¬ 
cus in Fig. Moreover, the absence of obvious eclipses, 
coupled with significant ellipsoidal-type photometric vari¬ 
ability, suggests 60° > i > 45°. Finally, it seems reasonable 
to suppose J-ijTx < I. The areas marked in grey in Fig. 13 
meet these three constraints. 

For each system, solutions that move the secondary into 
this grey zone can be achieved by increasing T 2 lT\\ by 
increasing both the mass ratio and reff( 2 ) over the default 
Ml values; or some combination of these. 

Implausibly large increases in secondary temperature 
are required to make the secondaries sufficiently luminous. 
Tolerable solutions are, however, possible by adopting val¬ 
ues of / T\ that are somewhat larger than those inferred 


This is true even shortly after leaving the main sequence, as 
an isolated massive star evolves to the right in the Hertzsprung- 
Russell diag ram at almost constant mass and luminosity (cf., 
e.g., Fig. pi. 
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Fig. 13: Constraints in the mass-luminosity plane for the primary and secondary components (upper, lower panels) of 
VFTS 450 and VFTS 652 (left, right panels), obtained by assuming that the primaries fill their Roche lobes, with absolute 
magnitudes and effective temperatures as summarized in Tablerefer to Section 5.5.1| for further details. 

Thick black solid lines show M-L loci for the indicated mass ratios q over a range in continuum brightness ratio, TijTx 
(values shown in the vertical scales to the right in each panel, marked at steps of 0 . 1 ; log(L/L 0 ) is constant for given 
TijJ-i, for fixed values). Thin solid and dashed curves, labelled Z and T in the top-lef t panel, show the zero-age 
and terminal-age main sequence loci for non-rotating single stars (from Brott et al.||2011a|b ). Thin red curves are lines 
of constant inclination, at i = 90, 60, and 45° (left to right). Grey shaded areas in the lower panels indicate the zones 
for which 45° < i < 60°, L > L(ZAMS), and TijTx < 1. 

Filled circles show the initial parameter estimates summarised in Table (columns headed ‘Ml’). ‘Error bars’ in the 
upper panels, and horizontal error bars in the lower panels, show the effects of changing Teff(l) by ±1 kK (this affects 
the inferred secondary mass, but not its luminosity, all else fixed). Vertical error bars in the lower panels show the effect 
of varying Teff(2) by ±2 kK (which has no effect on secondary mass). 

Open circles represent equivalent M-L solutions from V-band photometry. Green circles show the effects of (arbitrarily) 
adjusting to bring the secondary masses to 32 Mq, hence into the grey shaded zones in this plane. 

[Note that any changes to effective temperatures or absolute magnitudes also change the loci of constant g, so that only 
M and L can be inferred from this diagram for T^g or M{Iq) values that differ from the reference solution; e.g., the 
V-band solutions have the same g, i values as the /c-band solutions.] 
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Fig. 14: Constraints in the mass-luminosity plane for secondary components for fixed g, reff(l )i and absolute magnitudes, 
obtained by assuming that the primaries fill their Roche lobes (cp. Fig. 13); refer to Section [5.5.2 for details. 

Thick black solid lines show M-L loci for the indicated secondary temperatures (in kK) over a range in continuum 
brightness ratio, Tijl-i (values shown in the lower horizontal scales, marked at steps of 0.2). The secondary mass M is 
constant for given J^ 2 /-Fi(for fixed q values), as is the orbital inclination (upper horizontal scales, labelled in degrees). 
Thin solid and dash ed curves show the zero-age and terminal-age main sequence loci for non-rotating single stars (from 
Brott et al.||2011^ |. Grey shaded areas indicate the zones for which 45° < i < 60°, L > L(ZAMS), and T 2 lT\ < 1. 
Filled circles show the initial parameter estimates summarised in Table(columns headed ‘Ml’). Open circles represent 
equivalent results for F-band photometry. 


from the spectroscopy; this is not unreasonable given the 
errors, and the possibility that the secondary absorption¬ 
line spectra are ‘veiled’ by circumstellar material. 

Although an increase in J-ijTx suggests a compensating 
decrease in the secondary effective temperature (Table [^, 
we don’t attempt to refine the details further, given the 
considerable uncertainties (and noting that reductions of 
only < 10% in Teg are permitted if, at i < 60°, the secon¬ 
daries are not to be underluminous for their masses), but 
merely conclude that consistency between observations and 
broad evolutionary considerations can be achieved by plau¬ 
sible adjustments to the flux ratios, while retaining values 
for other parameters that are close to those estimated spec¬ 
troscopically. 

5.5.2. Secondary temperatures, brightness ratios 

Secondary temperatures are particularly prone to uncer¬ 
tainty if the secondary spectra are not purely photospheric, 
so we also explicitly examine the consequences of treating 
Teff(2) as a variable. For each system, the basic inputs are 
again (i) the absolute magnitude; (ii) the orbital period; 
(iii) the primary’s orbital-velocity semi-amplitude; with ad¬ 
ditionally (iv) the mass ratio and (v) the primary temper¬ 
ature fixed at selected values. 

System characteristics are now fully specified by Teg{2) 
and TijJ-x (assuming that the primary fills its Roche lobe). 
The consequences of varying these parameters are illus¬ 


trated in Fig. Once again, the simplest way to migrate 
the secondaries into the ‘zone of plausibility’ (shown in grey 
in the Figure) is to increase T^jTx, and/or to increase 
Teff(2). We conclude that, most probably, the secondary 
components are brighter than is superhcially suggested by 
the spectra, but that their spectra are veiled by accretion 
disks. The likelihood is that they are then also larger, and 
closer to filling their Roche lobes, than suggested by the Ml 
parameters, which could reconcile system properties with 
the photometry. 


6. Summary 

We have presented new spectroscopy of the massive blue 
binaries VFTS 450 and VFTS 652. Well-determined or¬ 
bits are established for the spectroscopically more conspic¬ 
uous components in both systems (the ‘primaries’ in our 
notation; Table [^; we argue that these are the less mas¬ 
sive components, and that they fill their Roche lobes, with 
near-synchronous rotation. Model-atmosphere analyses of 
the primaries yield reasonably robust results (Table [^, 
demonstrating signihcant surface-nitrogen abundances in 
each case. 

The secondary spectra have been detected, although 
the inferred characteristics are considerably less well estab¬ 
lished. If these secondary spectra reliably reflect the pho¬ 
tospheric properties of the secondary stars, then they are 
associated with the hotter components. However, quanti- 
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tative models of the systems built on the spectroscopic 
results (Table have inconsistencies with the photo¬ 
metry, and with evolutionary considerations, as discussed in 
Sections [5(4| and [5(5| We suggest that the secondary spectra 
are contaminated by, or arise in, accretion disks, and have 
explored the consequences of relaxing the allowe d values for 
relevant ‘observed’ secondary parameters (Figs.l^ 14). 
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Table 2: Log of spectroscopic observations of VFTS 450. Orbital phases are with respect to the circular-orbit ephemeris 
in Table The velocity measurements are described in Section 3.1 


Setting 

MJD 

Orbital 

Primary radial velocities (km s 


Secondary 



Phase 

He I 

He I 

Sim 


RV (km s-i) 




A4387 

A4471 

A4552 

Mean 

s.e. 

He II A4200 

UVES 

54761.2275 

0.99 

466.9 

460.7 

468.6 

465.4 

2.4 


UVES 

54761.2491 

1.00 

483.1 

466.0 

474.6 

474.6 

5.0 


UVES 

54761.2772 

0.00 

468.2 

476.6 

455.1 

466.7 

6.3 


UVES 

54761.2995 

0.00 

469.6 

460.7 

457.7 

462.7 

3.6 


UVES 

54767.2738 

0.87 

401.9 

385.0 

413.6 

400.1 

8.3 


UVES 

54767.2954 

0.87 

381.6 

383.6 

399.3 

388.2 

5.6 


UVES 

54845.1671 

0.17 

347.7 

355.7 

324.0 

342.5 

9.5 


UVES 

54845.1888 

0.18 

349.1 

358.4 

319.8 

342.4 

11.6 


UVES 

54876.1214 

0.66 

131.1 

123.3 

117.6 

124.0 

3.9 


UVES 

54876.1430 

0.67 

128.4 

128.6 

116.3 

124.4 

4.1 


UVES 

55172.2807 

0.63 

116.2 

115.3 

115.0 

115.5 

0.4 


UVES 

55172.3024 

0.64 

118.9 

111.3 

99.4 

109.9 

5.7 


UVES 

55173.2985 

0.78 

288.2 

280.0 

283.8 

284.0 

2.3 


UVES 

55173.3201 

0.78 

293.6 

285.3 

287.7 

288.9 

2.5 


UVES 

55178.1555 

0.49 

53.9 

52.9 

49.7 

52.2 

1.3 


UVES 

55178.1772 

0.49 

52.5 

55.5 

57.4 

55.2 

1.4 


LR02 

54748.2657 

0.11 

399.2 

403.6 

387.6 

396.8 

4.8 I 


LR02 

54748.2873 

0.12 

404.6 

399.6 

386.3 

396.8 

5.5 

>134 

LR02 

54748.3115 

0.12 

393.8 

395.6 

377.2 

388.9 

5.8 

LR02 

54748.3332 

0.12 

410.0 

398.2 

381.1 

396.5 

8.4 , 


LR02 

54749.2121 

0.25 

246.2 


228.0 

237.1 

9.1 


LR02 

54749.2337 

0.25 

240.8 

236.2 

224.0 

233.7 

5.0 


LR02 

54837.1223 

0.01 

462.8 

466.0 

460.3 

463.0 

1.6 'I 

>157 

LR02 

54837.1440 

0.01 

468.2 

463.3 

462.9 

464.8 

1.7 J 

LR02 

54868.0461 

0.49 

43.0 

51.5 

28.0 

40.9 

6.9 ” 

>400 

LR02 

54868.0677 

0.50 

45.8 

47.5 

28.0 

40.4 

6.2 J 

LR02 

55112.3016 

0.93 

450.6 

443.4 

456.4 

450.2 

3.8 ” 

>175 

LR02 

55112.3232 

0.93 

454.7 

443.4 

456.4 

451.5 

4.1 

LR02 

56210.3545 

0.25 

243.6 

236.4 

223.7 

234.6 

5.8 


LR02 

56210.3665 

0.25 

243.4 

236.4 

214.6 

231.5 

8.7 


LR02 

56210.3785 

0.25 

245.5 

230.5 

230.1 

235.4 

5.1 


LR02 

56217.3299 

0.26 

219.8 

211.2 

212.1 

214.4 

2.7 


LR02 

56217.3419 

0.26 

221.7 

207.5 

198.1 

209.1 

6.9 


LR02 

56217.3538 

0.26 

217.3 

206.2 

219.5 

214.3 

4.1 


LR02 

56243.3377 

0.03 

456.1 

457.2 

435.2 

449.5 

7.2 1 


LR02 

56243.3497 

0.03 

456.0 

458.8 

461.5 

458.8 

1.6 

U40 

LR02 

56243.3616 

0.04 

455.5 

457.8 

438.0 

450.4 

6.3 


LR02 

56256.2607 

0.91 

431.4 

415.0 

445.6 

430.7 

8.8 ’ 


LR02 

56256.2727 

0.91 

432.6 

414.1 

436.0 

427.6 

6.8 

lll5 

LR02 

56256.2846 

0.91 

425.5 

412.7 

440.4 

426.2 

8.0 


LR02 

56257.1301 

0.03 

460.9 

463.2 

457.9 

460.7 

1.5 ’ 


LR02 

56257.1421 

0.03 

465.9 

460.0 

469.2 

465.0 

2.7 

[ 86 

LR02 

56257.1541 

0.04 

463.1 

460.8 

430.1 

451.3 

10.6 


LR02 

56277.3081 

0.96 

474.3 

459.5 

487.0 

473.6 

7.9 ’ 


LR02 

56277.3201 

0.96 

467.9 

455.6 

469.3 

464.3 

4.4 

ll46 

LR02 

56277.3320 

0.96 

475.2 

459.4 

493.2 

475.9 

9.8 


LR02 

56283.0488 

0.79 

301.1 

288.0 

293.2 

294.1 

3.8 


LR02 

56283.0608 

0.80 

304.0 

284.3 

297.5 

295.3 

5.8 


LR02 

56283.0728 

0.80 

305.5 

287.2 

302.6 

298.4 

5.7 


LR02 

56294.1990 

0.41 

90.8 

79.8 

99.7 

90.1 

5.8 I 


LR02 

56294.2134 

0.41 

87.1 

82.7 

106.6 

92.1 

7.3 

I392 

LR02 

56294.2254 

0.42 

88.1 

81.9 

94.1 

88.0 

3.5 


LR02 

56295.1816 

0.55 

63.5 

63.5 

61.0 

62.7 

0.8 ’ 


LR02 

56295.1935 

0.56 

67.1 

64.0 

61.9 

64.3 

1.5 

U37 

LR02 

56295.2055 

0.56 

68.9 

67.9 

44.8 

60.5 

7.9 


LR02 

56304.2360 

0.87 

378.0 

366.7 

381.7 

375.5 

4.5 


LR02 

56305.2315 

0.01 

466.2 

461.2 

428.2 

451.9 

11.9 I 


LR02 

56305.2435 

0.01 

473.3 

464.0 

468.8 

468.7 

2.7 

h73 

LR02 

56305.2555 

0.02 

472.3 

459.7 


466.0 

3.6 

1 
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Table 2: continued 


Setting 

MJD 

Orbital 

Phase 

Primary radial velocities (km s 
He I He I Sim 

A4387 A4471 A4552 Mean 

s.e. 

Secondary 

RV 

(kms-i) 

LR02 

56306.2189 

0.16 

367.7 

366.3 

359.7 

364.6 

2.5 

LR02 

56306.2309 

0.16 

350.7 

367.7 

317.3 

345.2 

14.8 

LR02 

56306.2429 

0.16 

356.9 

360.8 

339.8 

352.5 

6.4 

LR02 

56308.1546 

0.44 

80.2 

69.0 

119.6 

89.6 

15.3 I 

1 

LR02 

56308.1666 

0.44 

87.2 

71.3 

87.7 

82.1 

5.4 1 

b23 

LR02 

56308.1786 

0.44 

84.0 

66.9 

71.8 

74.2 

5.1 

1 

LR02 

56316.2052 

0.60 

94.1 

87.6 

76.0 

85.9 

5.3 ’ 


LR02 

56316.2172 

0.61 

92.4 

92.6 

81.1 

88.7 

3.8 

b32 

LR02 

56316.2292 

0.61 

94.0 

91.6 

79.0 

88.2 

4.7 

I 

LR02 

56347.0132 

0.07 

444.9 

441.1 

420.1 

435.4 

7.7 ’ 


LR02 

56347.0251 

0.08 

441.2 

438.3 

462.3 

447.3 

7.6 

h32 

LR02 

56347.0371 

0.08 

437.6 

434.3 

416.1 

429.3 

6.7 

I 

LR02 

56349.0214 

0.37 

109.2 

88.5 

105.1 

100.9 

6.3 ’ 


LR02 

56349.0334 

0.37 

107.5 

92.8 

108.1 

102.8 

5.0 

b47 

LR02 

56349.0453 

0.37 

99.8 

91.9 

95.1 

95.6 

2.3 1 

I 

LR02 

56352.0241 

0.80 

289.7 

271.1 

279.1 

280.0 

5.4 ^ 


LR02 

56352.0360 

0.80 

291.9 

269.7 

282.6 

281.4 

6.4 

LR02 

56352.0480 

0.80 

296.8 

278.4 

280.4 

285.2 

5.8 

LR02 

56356.0044 

0.38 

90.7 

83.9 

91.5 

88.7 

2.4 I 

1 

LR02 

56356.0163 

0.38 

95.3 

81.2 

90.8 

89.1 

4.2 1 


LR02 

56356.0283 

0.38 

92.6 

79.0 

65.4 

79.0 

7.9 I 

1 




Si III 

He I 

He I 






A4552 

A4713 

A4922 



LR03 

54755.1987 

0.12 

387.0 

399.3 

402.7 

396.3 

4.8 

LR03 

54810.2265 

0.10 

412.4 

423.2 

433.8 

423.1 

6.2 

LR03 

54810.2481 

0.11 

410.5 

423.7 

425.3 

419.8 

4.7 

LR03 

54810.2699 

0.11 

409.6 

422.1 

419.7 

417.1 

3.8 

LR03 

54810.2915 

0.11 

404.9 

415.6 

419.6 

413.4 

4.4 

LR03 

54810.3248 

0.12 

395.2 

414.2 

412.6 

407.3 

6.1 

LR03 

54810.3465 

0.12 

394.2 

427.0 

413.5 

411.6 

9.5 
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Table 3: Log of spectroscopic observations of VFTS 652; details are as for Table 


Setting 

MJD 

Orbital 

Primary radial velocities (km s 


Secondary 



Phase 

He I 

He I 

Sim 


RV (km s-^) 




A4387 

A4471 

A4552 

Mean 

s.e. He II A4541 

UVES 

54791.3055 

0.69 

179.8 

171.1 

160.4 

170.4 

5.6 


UVES 

54791.3271 

0.69 

173.0 

172.4 

179.9 

175.1 

2.4 


UVES 

54792.1644 

0.79 

293.6 

299.9 

316.2 

303.2 

6.7 


UVES 

54792.1860 

0.79 

301.7 

293.3 

320.1 

305.0 

7.9 


UVES 

54847.1501 

0.19 

335.5 

321.2 

331.8 

329.5 

4.3 


UVES 

54847.1717 

0.19 

322.0 

317.2 

320.1 

319.8 

1.4 


UVES 

54892.0909 

0.42 

93.2 

91.4 

67.4 

84.0 

8.3 


UVES 

54892.1125 

0.43 

113.5 

91.4 

82.6 

95.8 

9.2 


UVES 

54894.0276 

0.65 

139.2 

135.2 

142.3 

138.9 

2.0 


UVES 

54894.0492 

0.65 

137.8 

132.6 

125.4 

131.9 

3.6 


UVES 

54896.0284 

0.88 

399.2 

398.2 

398.0 

398.5 

0.4 


UVES 

54896.0504 

0.88 

397.8 

396.9 

399.3 

398.0 

0.7 


UVES 

54897.0255 

1.00 

457.4 

480.1 


468.7 

11.3 


UVES 

54898.0463 

0.12 

408.6 

408.9 

400.6 

406.0 

2.7 


UVES 

54898.0687 

0.12 

401.9 

404.9 

377.1 

394.6 

8.8 


UVES 

55201.1762 

0.41 

91.8 

87.4 

78.7 

85.9 

3.9 


UVES 

55201.1984 

0.41 

89.1 

88.7 

94.2 

90.7 

1.8 


LR02 

54804.0935 

0.18 

343.7 

342.4 

346.1 

344.1 

1.1 1 


LR02 

54804.1151 

0.18 

338.2 

331.8 

347.4 

339.1 

4.5 


LR02 

54804.1368 

0.18 

338.2 

339.8 

344.8 

340.9 

2.0 

>232 

LR02 

54804.1584 

0.19 

336.9 

335.8 

338.3 

337.0 

0.7 

LR02 

54804.1801 

0.19 

332.8 

333.1 

334.4 

333.4 

0.5 


LR02 

54804.2016 

0.19 

324.7 

334.5 

329.2 

329.4 

2.8 


LR02 

54836.2280 

0.92 

431.7 

431.4 

440.8 

434.6 

3.1 1 


LR02 

54836.2497 

0.92 

429.0 

430.1 

438.2 

432.4 

2.9 

>212 

LR02 

54836.2758 

0.92 

430.3 

430.1 

439.5 

433.3 

3.1 

LR02 

54836.2974 

0.93 

442.5 

436.8 

436.9 

438.7 

1.9 . 


LR02 

54867.0978 

0.51 

60.6 

59.5 

52.7 

57.6 

2.5 ^ 

>329 

LR02 

54867.1195 

0.52 

63.4 

63.5 

61.8 

62.9 

0.5 

LR02 

55108.3179 

0.60 

83.7 

86.1 

82.6 

84.1 

1.0 


LR02 

55114.3094 

0.29 

183.9 

180.4 

186.4 

183.6 

1.7 'I 

>258 

LR02 

55114.3310 

0.30 

183.9 

179.0 

191.6 

184.8 

3.7 

LR02 

56210.3545 

0.90 

419.5 

410.1 

418.9 

416.2 

3.0 


LR02 

56210.3665 

0.90 

417.2 

404.2 

424.1 

415.2 

5.8 

I247 

LR02 

56210.3785 

0.90 

410.7 

411.4 

426.0 

416.0 

5.0 


LR02 

56217.3299 

0.71 

192.7 

182.5 

195.4 

190.2 

3.9 ’ 


LR02 

56217.3419 

0.71 

196.0 

181.8 

199.9 

192.6 

5.5 

I338 

LR02 

56217.3538 

0.71 

192.2 

185.4 

203.8 

193.8 

5.4 


LR02 

56243.3377 

0.74 

229.1 

212.2 

236.8 

226.0 

7.3 ’ 


LR02 

56243.3497 

0.74 

232.2 

215.3 

237.0 

228.2 

6.6 

I278 

LR02 

56243.3616 

0.74 

233.7 

217.6 

237.8 

229.7 

6.2 


LR02 

56256.2607 

0.24 

267.8 

240.1 

266.4 

258.1 

9.0 ’ 


LR02 

56256.2727 

0.24 

260.6 

236.7 

269.7 

255.7 

9.8 

[297 

LR02 

56256.2846 

0.24 

256.4 

242.9 

264.0 

254.4 

6.2 


LR02 

56257.1301 

0.34 

151.9 

128.2 

145.2 

141.8 

7.1 ’ 


LR02 

56257.1421 

0.34 

142.8 

131.9 

139.3 

138.0 

3.2 

I297 

LR02 

56257.1541 

0.34 

138.4 

128.7 

142.2 

136.4 

4.0 


LR02 

56277.3081 

0.69 

174.9 

173.2 

187.8 

178.6 

4.6 ’ 


LR02 

56277.3201 

0.69 

183.0 

172.7 

192.4 

182.7 

5.7 

I376 

LR02 

56277.3320 

0.69 

187.1 

177.6 

180.4 

181.7 

2.8 


LR02 

56283.0488 

0.36 

137.7 

133.8 

144.1 

138.5 

3.0 “ 


LR02 

56283.0608 

0.36 

137.3 

126.7 

129.7 

131.2 

3.2 

I348 

LR02 

56283.0728 

0.36 

135.2 

133.3 

132.2 

133.6 

0.9 


LR02 

56294.1990 

0.66 

145.7 

149.0 

152.5 

149.1 

2.0 ’ 


LR02 

56294.2134 

0.66 

149.4 

148.5 

150.0 

149.3 

0.4 

U94 

LR02 

56294.2254 

0.66 

153.9 

149.8 

147.3 

150.3 

1.9 


LR02 

56295.1816 

0.77 

278.2 

261.0 

291.5 

276.9 

8.8 ’ 


LR02 

56295.1935 

0.77 

277.7 

264.3 

283.6 

275.2 

5.7 

I245 

LR02 

56295.2055 

0.77 

283.2 

264.8 

294.7 

280.9 

8.7 


LR02 

56304.2360 

0.83 

342.8 

332.2 

357.8 

344.3 

7.4 

253 

LR02 

56305.2315 

0.94 

444.5 

438.9 

436.3 

439.9 

2.4 I 


LR02 

56305.2435 

0.94 

450.1 

442.5 

449.7 

447.4 

2.5 

U56 

LR02 

56305.2555 

0.94 

449.1 

442.7 

458.6 

450.1 

4.6 

1 
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Table 3: continued 


Setting 

MJD 

Orbital 

Primary radial velocities (km s 



Secondary 



Phase 

He I 

He I 

Sim 



RV (km s"^) 




A4387 

A4471 

A4552 

Mean 

s.e. 

He II A4541 

LR02 

56306.2189 

0.06 

449.4 

444.9 

457.9 

450.7 

3.8 


1 

LR02 

56306.2309 

0.06 

448.5 

449.5 

454.0 

450.7 

1.7 

1 

}215 

LR02 

56306.2429 

0.06 

451.8 

446.3 

453.1 

450.4 

2.1 


I 

LR02 

56308.1546 

0.28 

196.2 

175.1 

201.3 

190.9 

8.0 



LR02 

56308.1666 

0.28 

196.8 

178.2 

190.6 

188.5 

5.5 



LR02 

56308.1786 

0.28 

198.7 

178.4 

204.7 

193.9 

8.0 


I 

LR02 

56316.2052 

0.22 

300.9 

267.0 

293.7 

287.2 

10.3 

> 


LR02 

56316.2172 

0.22 

292.0 

268.7 

290.9 

283.9 

7.6 


hse 

LR02 

56316.2292 

0.22 

292.9 

269.1 

282.8 

281.6 

6.9 


I 

LR02 

56347.0132 

0.81 

319.8 

304.4 

327.5 

317.2 

6.8 



LR02 

56347.0251 

0.81 

316.2 

311.0 

325.4 

317.5 

4.2 


}255 

LR02 

56347.0371 

0.81 

321.1 

311.9 

326.4 

319.8 

4.2 


I 

LR02 

56349.0214 

0.04 

449.7 

447.3 

463.5 

453.5 

5.0 

> 


LR02 

56349.0334 

0.04 

450.8 

442.9 

459.8 

451.2 

4.9 


hoi 

LR02 

56349.0453 

0.04 

457.4 

448.2 

454.0 

453.2 

2.7 


I 

LR02 

56352.0241 

0.39 

105.5 

93.6 

111.1 

103.4 

5.2 

3 


LR02 

56352.0360 

0.39 

103.0 

99.4 

107.0 

103.1 

2.2 


b67 

LR02 

56352.0480 

0.39 

103.2 

99.7 

96.9 

99.9 

1.8 


I 

LR02 

56356.0044 

0.85 

370.5 

368.2 

379.4 

372.7 

3.4 

> 


LR02 

56356.0163 

0.85 

380.6 

367.4 

386.5 

378.2 

5.6 


}274 

LR02 

56356.0283 

0.86 

379.5 

368.2 

384.5 

377.4 

4.8 

J 

1 




Si III 

He I 

He I 








A4552 

A4713 

A4922 





LR03 

54808.1322 

0.65 

135.0 

136.2 

136.3 

135.8 

0.4 



LR03 

54808.1538 

0.65 

138.1 

140.6 

140.2 

139.6 

0.8 



LR03 

54808.1755 

0.65 

142.6 

143.2 

142.4 

142.7 

0.2 


>340 

LR03 

54808.1971 

0.66 

142.6 

149.8 

144.3 

145.6 

2.2 


LR03 

54808.2189 

0.66 

145.9 

146.6 

149.4 

147.3 

1.1 



LR03 

54808.2405 

0.66 

149.0 

150.9 

153.3 

151.1 

1.2 

_ / 
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